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We propose to probe the quantum ground state of a spin-1 Bose-Einstein condensate with the 
transmission spectra of an optical cavity. By choosing a circularly polarized cavity mode with an 
appropriate frequency, we can realize coupling between the cavity mode and the magnetization of the 
condensate. The cavity transmission spectra then contain information of the magnetization statistics 
of the condensate and thus can be used to distinguish the ferromagnetic and antiferromagnetic 
quantum ground states. This technique may also be useful for continuous observation of the spin 
dynamics of a spinor Bose-Einstein condensate. 

PACS numbers: 03.75.Mn, 32.10.Dk, 37.10.Vz, 42.50.Pq 



Spinor Bose-Einstein condensates (BECs) distinguish 
themselves from their scalar counterparts by their much 
richer internal structures and dynamics [l|, HJ. Based 
upon the single mode approximation (SMA), which is 
valid in most cases, the quantum eigenstates of a spinor 
BEC has been extensively and thoroughly studied with 
sophisticated algebraic techniques [3, 0, @]. Especially, 
for the spin-1 case, which we focus on in this paper, it 
is found that there are two possibilities for the ground 
state. Depending on a single parameter, the ground state 
can be either of an antiferromagnetic (AFM) or a ferro- 
magnetic (FM) nature. The two phases differ in their 
rotational properties — the AFM ground state is a SU (2) 
singlet and thus isotropic [||; while the FM ground state 
is a SU(2) multiplet heavily degenerate and thus can 
be well directed (up to the uncertainty principle associ- 
ated with the SU(2) algebra). From the point of view 
of quantum phase transition, these two phases deserve to 
be distinguished experimentally for their own sake. 

In this paper we propose a method based on cavity 
quantum electrodynamics techniques to fulfill this object. 
The basic idea is similar to that of Mekhov et al. in Ref. 

, where they proposed to probe the superfluidity- Mott 
insulator transition of cold atoms in optical lattices with 
cavity transmission spectra. In the dispersive regime, 
atoms inside the cavity shift the cavity frequency and in 
turn influence the cavity transmission spectra. The su- 
perfluidity and Mott-insulator phases differ in their atom 
number statistics and thus lead to drastically different 
cavity transmission spectra, which can be used inversely 
to differentiate the two phases. In their proposal, it is ad- 
equate to treat the atoms in a two-level fashion with the 
atomic levels structureless. However, here in our case, 
we have to go beyond the two-level model and take into 
account the detailed and realistic atomic level structure. 
In this way, we can retrieve the vector polarizability [8] 
of the atom that is dropped artificially in the two-level 
model. It follows that by choosing a circularly polarized 
cavity mode, the cavity photon couples not only to the 
total atom number but also to the total magnetization. 



The latter effect allows us to do our job, since the cavity 
transmission spectra carry information of the magnetiza- 
tion statistics of the BEC, which are absolutely different 
in the two phases. 

We assume a spin-1 BEC of N alkali atoms of a certain 
type is loaded inside an optical cavity [9j. This system 
is probed with a <r+ polarized laser so that the relevant 
cavity mode coupled with the atoms is also <r + polar- 
ized. The bare frequency to c of this cavity mode is cho- 
sen to resolve the alkali D line (fine structure) and placed 
relatively near the D\ transition [Toj ] . However, the de- 
tuning between the D\ transition and the cavity mode 
frequency is still large enough not to resolve the hyper- 
fine structures of the D\ line. More precisely, we have 
A ff Fsi < |w c - um\ < |w c - £Jd 2 |, where A H fsi is the 
magnitude of the hyperfine splitting of the D\ line, ujdi 
and u)D2 are the average frequency of the hyperfine lev- 
els of D\ and D2 , respectively. In the appendix we show 
that the effective Hamiltonian describing the atom-cavity 
photon interaction is (H = 1 throughout in this paper) 

Hint = ^ f 1 + y) Uo J ' dxu 2 {x)ii)l n {x)i}j m {x)a { a. (1) 
m 

Here ip m (x) (to = 0,±1) is the field annihilation oper- 
ator for an atom in the hyperfine state \F — 1, to) at 
position x (the quantization direction is along the cav- 
ity axis and is denoted as the z-dircction), and a is the 
annihilation operator for the cavity photon. The cavity 
mode function with its maximum normalized to unity is 
denoted as u(x). The parameter Uq characterizes the 
strength of the ac-Stark shifts of the atomic levels per 
photon, in the viewpoint of the atoms; or alternatively, 
the magnitude of the cavity frequency shift per atom, in 
the viewpoint of the photons. Formally, we can write Uq 
as Uq = <?q/ '(cj c — Wei) by introducing an effective atom- 
photon coupling strength go, which incorporates the over- 
all effect of coupling with all the hyperfine levels in the 
D\ line. 

Here some comments are in order. The main feature of 
Hamiltonian (fTJ) is the spin dependence of the level shifts, 



r 
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which can be interpreted in terms of a fictitious magnetic 
field. In most experiments aiming to study intrinsic spin 
dynamics of a spinor BEC, this effect is (and should be) 
avoided or minimized by taking the laser beam to be lin- 
early polarized or with a detuning too lar ge t o resolve the 
D line or both (an exception being Ref. However, 
for the purpose of continuous observation of spin dynam- 
ics and spin magnetization, this effect proves to be useful. 
Actually, it is exactly the same physics that underlies the 
Faraday rotation spectroscopy develo ped recently for in 
situ observation of a spinor BEC (l2l. Il3l. Il4j . For if the 
cavity mode is <r_ polarized, the plus sign in front of ^ 
in Eq. ([I]) should be replaced with a minus sign. That 
is, the cr+ and er_ polarized light respond differently to a 
given atomic sample — the essence of the Faraday rotation 
effect. 

For the spinor BEC, we assume the SMA is valid and 
rewrite the atomic field operators as ip m (x) = (f>(x)c m . 
Here 4>(x) denotes the common spatial mode (normal- 
ized to unity) for the three spin components and c m is the 
annihilation operator associated with the corresponding 
spin component. Under the SMA, we can rewrite Hamil- 
tonian as 



Hint = U N 
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where N = c ln c m is the total atom number and 
M = J2 m mc ln c m is the total magnetization. The pa- 
rameter U = UqJ dxu 2 (x)(j) 2 (x) is proportional to the 
overlap between the cavity mode function and the con- 
densate mode function. In general, the extension of a 
condensate is much larger than the period of u 2 (x) in the 
axial direction while much smaller than its waist in the 
transverse direction. Therefore, for a condensate placed 
around the axis of the cavity, the integral can be well 
approximated with 1/2, i.e., U ~ Uq/2. We note that 
N and M commute with both Hi nt and the Hamiltonian 
of the spinor BEC itself. Therefore, we have a quantum 
non-demolition measurement. 

Assume the probe laser is of frequency lu p . The Heisen- 
berg equation (in the frame rotating at lo p ) for the cavity 
mode annihilation operator is 

a = —i ^oj c — lo p + U (^N + 2"^N ^ a — Ka + r), (3) 

where k is the damping rate of the cavity mode and rj 
the driving amplitude. The damping term leads to the 
stationary solutions of a and the photon number Q, 
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where A = u> c + UN — ui p is the shifted cavity-probe 
detuning. The expectation value of n st over a prescribed 



state of the spinor BEC is then 
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where Jm is the probability for M taking value M in the 
prescribed state. We have J2ai = 1' wnere the sum 
is over all the possible values of M. As pointed out in 
Ref. [HI], the total photon number contains incoherent 
contributions of a scries of coherent states, the M-th of 
which is of amplitude i]/(K + i(A + UAl/2)) and of weight 
fnj. Mathematically, we see from Eq. ([5]) that n s t, as a 
function of the detuning A, is a weighted sum of a series 
of Lorentz functions, which are evenly displaced with U/2 
and of full-width half-maximum 2k. The point is that the 
cavity transmission spectra (proportional to n st ) contain 
information of the atomic statistics embodied in /^/. 

We now apply the formalism developed above to the 
two possible phases of a spin-1 BEC. For a spin-1 BEC 
in the SMA, the Hamiltonian accounting for the spin 
dynamics is H s = W'S 2 Q, with the spin operator 
S = (S x ,S y ,S z ) = 4(5™",5™",^ m ")c™- The param- 
eter W is proportional to the difference of the s-wave 
scattering lengthes in the two total spin channels (F = 
or 2). Obviously, the sign of W determines the nature of 
the ground state. If W > 0, the ground state is of total 
spin 5 = and given explicitly as 0, H, @ 
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= \S = 0,M = 0) oc (A^) N/2 \0), (6) 



where the scalar (rotationally invariant) operator A 1 
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This state is a singlet and also the unique 
ground state. For this AFM ground state, the magneti- 
zation has no fluctuations at all. The only nonvanishing 
term in Eq. ([5|) is M = 0. Therefore, by scanning the 
probe frequency u> pi we get the intracavity photon num- 
ber as a single Lorentz function, 
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with a width of 2k and maximum height of unity (in 
units of rj 2 / k 2 , the maximum possible intracavity photon 
number). Note that it is of the same shape as that of a 
bare cavity, though the center is shifted by UN. 

If W < 0, the ground state is of total spin S = N and 
2N + 1 fold degenerate. Explicitly, we have as a set of 
orthonormal basis [1, [f| , 

\G) fm = \S = N, M) ex (5_) Ar - M (4 1 ) A, |0), (8) 



where S- — v^cjc+i + c^co) is the lowering operator 
for S z . These 2N + 1 states each have a definite mag- 
netization, thus the corresponding transmission spectra 
is also a single Lorentz function (but the center of which 
may vary). However, here we are most interested in the 
experimentally most accessible states — angular momen- 
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turn coherent states — states that can be obtained by ro- 
tating the z-aligned state \N,N), 

e-^ § \N, N) = {X+1C ^ + X - lC - l)JV 10). (9) 



Here C — ( — sin (j), cos (j), 0) perpendicular to the z-axis 
is the axis of rotation, while 9 is the angle of rota- 
tion or inclination angle. The Xi S clXG determined by 
e -iec-s c ^ e iec-s _ a ; +lC t_ i + XoC t + x-xcLi- Tne state 
constructed in Eq. ([9]) defines a direction (9, <fi) in that it 
is the eigenstate of S n ^ s ^) = sin 9 cos <j>S x +sin sin + 
cos 9S Z with eigenvalue 2V. These states can be readily 
prepared with a proper radio-frequency (rf ) pulse from 
the state \N, N) For such a state, the weight coeffi- 
cients of the Fock states \N + i, Nq, N-i) follow the triple- 
nominal distribution 



PiN+uNo^) = ^^ pZtWfr, (10) 

with pi = \xi\ 2 and J2iPi = 1- Explicitly, we have 
(p + i,po,P-i) = (cos 4 | , 2 cos 2 | sin 2 | , sin 4 | ) . Note that 
Pi are independent of the azimuth angle <f>. This is rea- 
sonable since all angular momentum coherent states with 
the same inclination angle 9 are connected by rotations 
around the z-axis. Therefore they share the same distri- 
bution P(N + i,N ,N-i). 

The probability distribution function /m can be gen- 
erated in a random walk process. Image a particle, ini- 
tially at the origin, performing random walks on the 
real axis. At each step, the particle either moves right- 
ward or leftward or just keeps still, with probabilities 
P+it P-i, an d po, respectively. Then /m is just the 
probability of arriving at M after N steps. In other 
words, M is the sum of N independent random variables 
{Xk,l < k < N}, the common distribution function of 
which is Pr(Xk = i) — Pi, with i = 0, ±1. This in- 
terpretation enables us to calculate /m both analytically 
and numerically. Analytically, for N S> 1, by the central 
limit theorem, we can approximate Jm with the normal 
distribution 
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where M c — N cos 9 and a = y/~N do (here Co = y - 
being the variance of Xf.) are respectively the average 
value and variance of M. Numerically, the characteristic 
function of f M is x(*) = (p+ie +it + Po +P-ie~ it ) N . By 
an inverse Fourier transform, we get /m determined (it is 
just the coefficient of e iMt in the expansion of x(i)). The 
two approaches yield the same result up to a high preci- 
sion. In Fig. [TJa), we show the distribution Jm for four 
different inclination angles with the total atom number 
TV = 10 6 . We see the larger the absolute value of sin6>, 
the larger the fluctuations of M. For all cases, we see a 
clear y^/V scaling of a . 



We calculate the intracavity photon number n st by ap- 
proximating /a/ with Eq. pip and replacing the sum with 
an integral. The integral can be further understood as 
the convolution of the probability distribution function 
Jm and a Lorentz function. We can then use the con- 
volution theorem to first solve the Fourier transforms of 
the two functions respectively, form their product, and 
then take an inverse Fourier transform. We get 



n st 
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(A) = IL / dte u( - A ~ A ^ x e- K l ; 

2 K J-oo 



',(12) 



where A c = —M c U/2. Obviously, the procedure can be 
done inversely, i.e., from the cavity transmission spectra 
observed experimentally, we can infer the distribution of 
the magnetization by Fourier transforms. 

We can read two facts from Eq. (fT2"]) . Firstly, the shape 
of n st as a function of A is shifted by A c from that de- 
picted in Eq. ([7]). This overall shift is due to the average 
deviation of M from zero and thus has a classical na- 
ture. Secondly, the shape of n st is broadened compared 
to that in Eq. ([7]), which is due to the fluctuations of M 
around the mean and thus has a quantum nature. The 
shift and broadening are appreciable if (M C U, aU) ^> k, 
or in terms of N, (NU, VNU) » k. Currently, g a is 
about 2ir x 15 MHz and the line width k of the cavity 
can be made as small as 2-k X 1.0 MHz [13, E El- By 
taking the detuning A ca — uj c — lujji to 2tt x 40 GHz, 
we get U ~ 2tt x 2.8 kHz. The atom number of a spinor 
BEC is typically on the order of 10 6 . Thus we see the 
first condition NU ^> n is well satisfied while the second 
condition \/NU ^ k is marginally satisfied. In the near 
future, we expect re can be reduced by one order and then 
the second condition can be well met also. 

In Fig.[TJb), we present the intracavity photon number 
at steady state n st as a function of the detuning A. For 
the AFM ground state, the curve (solid line) is a sim- 
ple Lorentzian centered at A c = and with width 2re. 
This is the only possibility. However, for the FM ground 
states, depending on the inclination angle 9, the cen- 
ter and width of the curve vary widely (the figure may 
be a bit misleading, but note that we have artificially 
shifted the curves to a common center; see the labels of 
the horizontal axes). Especially, the broadening effect 
is apparent compared to the AFM case (here we have 
VNU/k = 10). Thus the AFM and FM ground states 
can be distinguished on one hand by the different shapes 
of the corresponding spectra [20j and on the other hand 
by their responses to rf pulses (or a transverse magnetic 
field). Under the action of rf pulses, the AFM ground 
state remains unaltered since it is rotationally invariant, 
so the cavity transmission spectra are unchanged. In 
contrast, a FM ground state can be tickled to other di- 
rections and thus both the center and the width of the 
transmission spectra will change. 

An alternative and more convenient method may be 
just reversing the polarization of the probe laser. Be- 
cause the (7+ and <r_ polarized cavity modes couple re- 
spectively to M and —M, it is easy to see that for the 
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FIG. 1: (Color online) (a) Probability distribution function 
/m of the magnetization M of ferromagnetic (FM) ground 
states with different inclination angles, (b) Intracavity photon 
number at steady state n at as a function of the detuning A. 
The solid line corresponds to the antiferromagnetic (AFM) 
ground state. Other lines correspond to FM ground states in 
(a). The parameters are N = 10 6 and U = 0.01k. For the 
horizontal labels, M c = NcoaO and A c = —M c U/2 in the 
FM case, while M c = and A c = in the AFM case. 



AFM case, the transmission spectra are unaffected; while 
for the FM case, the transmission spectra are reflected to 
be centered at — A c . The distance between the two cen- 
ters is on the order of NU, which is much larger than k, 
therefore this effect can be readily observed. This method 
makes evident the similarity between our approach and 
the Faraday rotation spectroscopy. 

In summary, we have shown that by choosing a cir- 
cularly polarized cavity mode, we can achieve coupling 
between the photon number with the magnetization, i.e. 
the z-component of the total spin of the spinor BEC. Full 
statistics of the magnetization of the spinor BEC can be 
extracted from the cavity transmission spectra and this 
allows us to discern the AFM and FM phases, which 
due to their different rotational properties, possess abso- 
lutely different magnetization statistics. Our discussion 
has been focused on the thermodynamical properties of 
the BEC (ground state oriented). However, it is clear 
that the technique can also be used for observing the 
real-time evolution of a spinor BEC. Moreover, since the 
condition NU ^> k can be readily fulfilled, the sensitivity 
can be high. 

It should be emphasized that our discussion has been 
greatly simplified by neglecting the center-of-mass mo- 
tion of the BEC and by the SMA [without the two 
approximations, we cannot reduce Eq. |T]) to Eq. ([2"])]. 
These approximations are reasonable if the probe is weak 



enough. For future works, however, the strong probe case 
is worth investigation. If the probe is so strong that the 
intracavity optical lattice cannot be neglected, center-of- 
mass motion of the spinor BEC can be excited. Moreover, 
since the optical lattice felt by different spin components 
are different, a dynamical SMA is inappropriate. The sit- 
uation is further complicated by the back-action of the 
BEC on the cavity mode. In the end, this system may 
exhibit many interesting nonlinear effects [2l| . 
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2006CB921206, and 2006AA06Z104. J. M. Z. would like 
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APPENDIX: EFFECTIVE HAMILTONIAN 

In this appendix we review the derivation of the ef- 
fective Hamiltonian ([T]). The main concern is to demon- 
strate the Zeeman-like effect. General theor y ab out this 
effect has actually been well established [22l. l23l |24T | . We 
include the derivation here just to serve our concrete pur- 
pose. For simplicity, we take the optical field to be clas- 
sical, generalization to the quantum case is straightfor- 
ward. The alkali atom in question is assumed to be in 
the n 2 Si/ 2 F = 1 hyperfine ground state manifold. A 

laser with electric field E{t) — E exp(—iujft) + c.c. cou- 
ples the atom to the n 2 P 1 / 2 excited state manifold (Di 
line). Here it is assumed that the frequency of the laser 
is chosen so that the fine structure of the alkali D line is 
well resolved while the hyperfine structure is not. This 
is feasible since the fine structure splitting is generally 
three or four orders larger in magnitude than the hyper- 
fine structure splitting for the alkali D line. For instance, 
the fine splitting of the D line of 87 Rb is 2tt x 7200 GHz, 
while the hyperfine splitting of the D\ line is 27r x 814.5 
MHz. 

We aim to derive an effective Hamiltonian confined in 
the F = 1 ground state subspace. By second order per- 
turbation theory and taking the rotating wave approxi- 
mation, we have 12511 



H, 



eff 



Pg;F=l{E* ■ d)P D1 {d- E)P g , F = 1 
LUf — UJ]Jl 



(A.l) 



where d is the electric dipole operator and P g -F=i 
and Pdi are the projection operators on the F = 
1 hyperfine ground state manifold and Dl fine ex- 
cited state manifold respectively. That is, P 9 -f=i = 
E„> 2 Si/25^ = l,m)(n 2 S 1/2 ;F = l,m\ and ' P D1 = 
J2Fm \ n2 Pi/2', F,m)(n 2 P 1 /2; F,m\, where the sum is over 
all the possible values of F and m. For the latter, we can 
also transfer from the coupled representation to the un- 
coupled representation, i.e., equivalently we have 



P, 



Dl 



(A.2) 
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Here P e ;J=i/2 is purely an electronic projection opera- 
tor onto the J = 1/2 excited state subspace while Pj is 
purely a nuclear projection operator. We can also define 
the project operator 



which projects onto the ground state manifold up to the 
fine structure. Substituting (|A.2|) and (|A.3p into (|A.1|) 
and noting that P g P g .p = i = P g -F=i, we rewrite the ef- 
fective Hamiltonian as 



Pn 



P g;J=l/2 



Pj 



Pg-F=l + Pg;F=2, 



(A.3) 
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P g (E* ■ d)P D1 (d ■ E)P g 



Pn 



g;F=i 



p ( P 9 ;J=iME* ■ d)P ei j= 1/2 (d ■ E)P g , J=1/2 \ 

Pg;F=l ® Pi Pg;F=\- 

\ Wf-w D1 J 

I 



(A.4) 



Note that in the parentheses we have successfully sepa- 
rated the electronic part from the nuclear part. This is 
important since the external electric field couples with 
the electrons but not the nucleus. We define D = 
Pe-,j=i/2dP g ;j=i/2i which is a vector operator (and so is 
D^). The electronic part in the parentheses in (|A.4|) is 



then in the form (E* ■ D^){D ■ E), which can be rewrit- 
ten in the form D\DjE* Ej, the scalar contraction of two 
rank-2 tensors, with the atomic and optical parts sepa- 
rated. We can decompose the dyadic DjDj into rank-0, 
1, 2 irreducible tensors, 
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with E*Ej done similarly. We note that since D'D 
acts on a two-dimensional Hilbert space, the rank-2 irre- 
ducible part vanishes identically. We can then contract 
the corresponding irreducible tensors to form scalars. 
The result is 



1 



& -D) E* -E 
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Z3 f x 3) ■ [E* x E) . (A.6) 



Using the Wigner-Eckart theorem, we see that in the 
J = 1/2 ground state subspace the scalar operator 3* ■ D 
is proportional to the identity operator while the vector 
operator D< x 3 to J. Thus the above equation is equiv- 
alent to 



\E\ 2 [a - i/3(e* x e) • fj 



(A.7) 



with a and (3 being some constants and e — E/\E\ being 
the laser polarization vector. Substituting this result into 



(fA~4ll and using again the Wigner-Eckart theorem, we 
have finally 
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eff 



LUf - LUD 
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a — i^(e* x e) • F ) , 



(A.8) 



with 7 being some constant incorporating the Lande fac- 
tor [261 ]. The first term describes a center-of-mass shift 
of the ground state multiplet while the second term is 
in effect a Zeeman term — the very effect we are seeking. 
Note that for a linearly polarized laser the second term 
vanishes identically while for a + and ct_ polarized lasers, 
the fictitious magnetic fields point opposite. 

Therefore we have shown how an effective Hamiltonian 
like (1) can arise naturally under certain ciucumstances. 
To get the coefficients a and 7 determined, it is bet- 
ter to resort to the dipole matrix elements tabulated in 
Ref. [13. 
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